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Abstract 

Some properties of cosmological models with matter creation are inves- 
tigated in the framework of the Friedman-Robertson- Walker (FRW) 
line element. For adiabatic matter creation, as developed by Prigogine 
and coworkers, we derive a simple expression relating the particle num- 
ber density n and energy density p which holds regardless of the mat- 
ter creation rate. The conditions to generate inflation are discussed 
and by considering the natural phenomenological matter creation rate 
tjj = 3f3nH, where /3 is a pure number of the order of unity and H 
is the Hubble parameter, a minimally modified hot big-bang model 
is proposed. The dynamic properties of such models can be deduced 
from the standard ones simply by replacing the adiabatic index 7 of 
the equation of state by an effective parameter 7^, = 7(1 — (3). The 
thermodynamic behavior is determined and it is also shown that ages 
large enough to agree with observations are obtained even given the 
high values of H suggested by recent measurements. 



1 Introduction 



The origin of the material content (matter plus radiation) filling the 
presently observed universe remains one of the most fascinating unsolved 
mysteries in cosmology even though many authors worked out to understand 
the matter creation process and its effects on the evolution of the universe 
[1-27]. 

Radiation and matter constituents can quantum-mechanically be pro- 
duced in the context of Einstein's theory or, more generally, in any relativistic 
theory of gravitation. Such a process has been systematically investigated by 
Parker and coworkers by considering the Bugoliubov mode- mixing tech- 
nique in quantum field theory. This approach, roughly speaking, follows 
naturally from the fact that in curved spacetimes, as well as in accelerated 
frames, it is usually impossible to fix a priori a unique vacuum state for 
quantum fields]^]. In particular, this means that an observer with a detector 
will detect at late times a nonvanishing fiux of particles in a state initially 
set up to be empty of particles (Fulling- Unruh-Davies effect) . Unfortunately, 
due to expected back-reaction effects, it is not so clear that such a mecha- 
nism can account for sufficient particle creation to explain either the cosmic 
background radiation or the matter content of the universe. 

An alternative approach to matter creation, the development of which 
gave rise to deep physical insights for theoretical cosmology, was suggested 
by TryonQ and independently by Fomin[Q. They argued that if the net 
value of all conserved quantities of the universe is zero, as for instance, the 
total energy (gravitational plus material), then a universe whose duration is 
quantically restricted by the uncertainty relation AEAt > h/2, could have 
emerged as a vacuum fluctuation. No speciflc scenario was proposed by these 
authors however, in such an approach, the universe should be spatially closed 
in order to have all net charges identically zero. These ideas guided Zeldovich 
and other researchers to investigate the possibility that the classical space- 
time came into existence from a quantum-gravitational tunneling process 
termed "spontaneous birth of the universe" [§, |[ . 

A different but somewhat related line of development was pursued by 
Brout, Englert and Gunzig|^. They proposed a concrete scenario that pro- 
vided a simultaneous generation of matter and curvature from a quantum 
fluctuation of the Minkowski spacetime vacuum. In this model, after a flrst 
stage of creation, the universe enters a de Sitter phase by which some cosmo- 
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logical problems (horizon, flatness, etc.) are solved. In a subsequent period, 
the system finally achieves the standard FRW phase. 

Many attempts to treat the matter creation process at a phenomenolog- 
ical macroscopic level have also long been considered in the literature based 
on rather disparate motivations (for a review of early literature see Ref . p|] ) . 
There have also been some claims [9-11], that particle creation during or near 
the Planck era could classically be modeled by bulk viscosity stress(second 
viscosity). This is an interesting connection since irreversible processes are 



believed to play a fundamental role in the problem of time-asymmetry [|12| . In 



this case, the usual thermodynamic "arrow of time" , translated in this con- 
text as entropy generation due to matter creation, could provide a natural 
explanation of the "arrow of time" in the cosmological domain. The enlarge- 
ment of the traditional FRW equilibrium equations to include these effects 
has also, at least, two additional goals, namely: to explain the observed high 
entropy of the cosmic background radiation and to avoid the initial singular 
state existing in the standard equations [13-16]. 

More recently, irreversible processes have become the subject of study 
once again in connection with infiationary universe scenarios [17-20]. The 
basic idea is that bulk viscosity (matter creation) contributes at the level of 
the Einstein field equations (EFE) as a negative pressure term. It turns out 
that effective negative pressures are the key condition to generate inflation. 



Speciflcally, Barrow [[181 introduced this idea in the framework of the new 
inflationary scenario. He claimed that particle creation due to nonadiabatic 
decay of the fleld driving the slow-rollover inflation can macroscopically be 
described by the viscous cosmological model found by Murphy [^, for which 
the bulk viscosity coefficient is proportional to the energy density of the fiuid. 

The above considerations show that since the very beginning bulk vis- 
cosity has been widely interpreted as a phenomeno logical description of the 
matter creation process in the cosmic fiuid(see Refs. quoted in [21] for recent 
papers in this line). However, regardless of this macroscopic analogy as well 
as any microscopic description, it is important in itself to know how matter 
creation can be incorporated in the classical Einstein field equations. This 
question was seriously considered in the pioneering article of Prigogine and 
coworkers |2^, who implicitly pointed out that the bulk viscosity and matter 
creation are not only independent processes but, in general, lead to different 
histories of the universe evolution. They argued that, at the expense of the 
gravitational field, matter creation can occur only as an irreversible process 
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constrained by the usual requirements of nonequilibrium thermodynamics. 
The crucial ingredient of the new approach is the explicit use of a balance 
equation for the number of created particles in addition to the Einstein field 
equations. When properly combined with the thermodynamic second law, 
such an equation leads naturally to a reinterpretation of the stress tensor 
corresponding to an additional negative pressure term which, as should be 
expected, depends on the matter creation rate. This is in marked contrast to 
the bulk viscosity formulation, in which entropy is produced but, the num- 
ber particle conservation law is taken for granted. These results were further 
discussed and generalized by Calvao, Lima and Wagap3|, through a co- 
variant formulation allowing specific entropy variation as usually expected 
for nonequilibrium processes in fluids. The issue of why the processes of bulk 
viscosity and matter creation are not equivalent either from a dynamic or a 
thermodynamic point of view has been recently discussed in the literature 

mm- 

The macroscopic irreversible approach to matter creation has also been 
applied by these authors to early universe physics. For instance, Prigogine et 
al[^ obtained a scenario quite similar to that proposed in Ref.0, in which 
the universe emerges from an initial Minkowski vacuum. Thus, to a cer- 
tain extent, this work can be viewed as the macroscopic counterpart of the 
ideas originally proposed by Tryon0] and Fomin[^] and further semi classi- 
cally developed in the model proposed by Brout et al.|0. However, unlike the 
reversible semiclassical equations considered in the latter, the phenomeno- 
logical approach provides, in a natural way, the entropy burst accompanying 
the production of matter. 

In this article we focus our attention on the "adiabatic" matter creation as 



originally formulated in Ref. and somewhat clarified in As we shall 



see, unlike the standard model, to construct a definite scenario with matter 
creation one needs to solve a system of three coupled differential equations 
since the balance equation for the number density has been added to the pair 
of independent EFE. In principle, the full integration of such a system is not 
a trivial task because it depends on the somewhat unknown matter creation 
rate. However, as will be seen, it is possible to establish a simple relation 
between the equations for the particle number and energy density which holds 
regardless of the matter creation rate. This result will allow us to write the 
differential equation for the scale factor in terms only of the matter creation 
rate, thereby simplifying the analysis of the physically admissible models as 
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well as their comparison with the bulk viscous universes. The conditions 
to generate inflation in this context will be generically discussed, however, 
unlike previous work on this subject [22-27], we are more interested in the 
late stages of universe evolution. In this sense, a new class of cosmologies 
endowed with matter creation, leading to definite predictions in the present 
phase, is proposed. As argued, for all values of the curvature parameter, 
this is the simplest class of hot big-bang cosmologies driven by the matter 
creation process. As in the standard model, the thermodynamic behavior 
is readily computed and it is also shown that ages large enough to agree 
with observation can be obtained even given the high values of the Hubble 



parameter suggested by the recent measurements |28, 29 



2 FRW Equations With Matter Creation 

Let us now consider the FRW line element (c = 1) 

ds^ = de - R\t){^^ + rHe^ + rhin\d)d(t>'') , (1) 
1 — kr^ 

where A; = 0, ±1 is the curvature parameter of the spatial section and R is 
the scale factor. 

In that background, the nontrivial EFE for a fluid endowed with matter 
creation and the balance equation for the particle number density can be 
written as [23-25] 

SttGp = 3^^ + 3-^ , (2) 
8.G(p + p.)^-2|-|-A . (3) 

^ + 3§ = ^ , (4) 

n R n 

where an overdot means time derivative and p, p, n and if) are the energy 
density, thermostatic pressure, particle number density and matter creation 
rate, respectively. The creation pressure Pc depends on the matter creation 
rate, thereby coupling Eqs. (3) and (4) to each other and, although indirectly. 
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both of them with (2) through of the energy conservation law which is con- 
tained in the EFE themselves. For "adiabatic" matter creation, this pressure 
assumes the following form (See Ref. EM for a more general expression) 



where H = R/R is the Hubble parameter. 

As it stands, the above system is underdetermined since there are six 
unknowns, namely: p, p, Pc, n, R, ip and only three equations plus the 
constraint (5). It thus follows that one needs to provide two more relations 
in order to construct a definite cosmological scenario with matter creation. 
The first constraint takes the form of an equation of state which is supplied 
by thermodynamical considerations. The one usually employed in cosmology 
is the so-called "gamma-law" equation of state 

P = (7 - 1)P , (6) 

where the constant 7 lies in the interval [0,2]. The second constraint is a 
specification of the exact form of the matter creation ip (which should be de- 
termined from a more fundamental theory involving quantum processes). At 
this point, the procedure followed in the literature has been: (a) to integrate 
Eq. (4) assuming a given phenomenological law for ip (b) to insert the expres- 
sion of n into (5) and using (6) to obtain the evolution equation for the scale 
factor. Here, we consider a more general and somewhat more comprehensive 
approach. Firstly, we will visualize the kind of coupling existing among the 
balance equation (4) and the EFE. To that end, we establish the differential 
equation for i? as a function of ip and n. Combining Eqs.(2) and (3) with 
(5) and (6) it follows that 

This is a very enlightening expression in determining the effects of ip on 
the evolution of the scale factor. Since the 3nH term in (4) measures the 
variation of n only due to the expansion of the universe, it proves convenient 
to introduce the dimensionless and in general time-dependent parameter 
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in order to measure the effects of tlie matter creation rate. As expected, for 
= P = 0, Eq.(7) reproduces tlie general FRW equation, tliereby decoupling 
the subsystem formed by equations (2) and (3) from (4). Indeed, unlike 
claims by the authors of [Q, such decoupling ior ip = happens only if the 
"7 - law" has been adopted (Ref. ||30| examined the couphng of (2) - (4) 



for = and a more general equation of state). If ijj is different from zero 
but ip « 3nH, that is, (3 « 1, the effects of ip may be safely neglected. 
Physically, one may expect that the most interesting solutions of (7) arise 
during the phases in which the parameter (3 is of the order of unity. 

To proceed further, we now establish a general expression relating n and 
p. By inserting (5) and using (6) in the energy conservation law 

p + 3ip + p + p,)H = , (9) 

it takes the form 

p + 3-fpH = -fp^ . (10) 
n 

Therefore, comparing (10) with (4) it follows that 



7- = - , (11) 

n p 



the solution of which is 



n = no{-y , (12) 

Po 

where rio and po are the values of n and p at a given instant (from now 
on the index o denotes the present values of the parameters). It is worth 
mentioning that (12) holds regardless of the specific form assumed for the 
matter creation rate tp. It reduces to the right limit for a dust filled universe 
since in this case p = nM, where M is the mass of the created particles. 
It is also easy to see that (12) does not remain valid in the more general 
formulation proposed in Ref. P^ . 

In summary, the set of equations (2)-(6) has been reduced to Eq.(7) to- 
gether with relation (12). Thus, taking into account relation (2) for the 
energy density p, the integration of the evolution equation for the scale fac- 
tor depends only on the form of ip. In principle, since matter creation is 
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essentially a quantum process, the corresponding rates should be obtained 
from a quantum field theory in the presence of gravitational fields. However, 
in the absence of a well accepted model, the natural way is to investigate 
physically interesting solutions of (7) by adopting a phenomenological de- 
scription. In the case of viscous models, for instance, the analogous step is to 
to prescribe a form of the bulk viscosity coefficient. The one widely adopted 
in the literature is ^ = ap'^ , where a is a dimensional constant and v lies in 
interval [0,1] (see Ref. 

For matter creation, Prigogine et al. [^] examined the consequences of 
assuming a rate ip = aH^ (a constant) for a dust filled FRW fiat model. As 
one can see from (12) and (2), such a choice corresponds to ip = aSiiGp/S, 
that is, the same type of phenomenological expression first considered by 
Murphy [|l^ in the context of a cosmology with viscosity. However, as dis- 
cussed in [^, the models are quite different from a physical point of view. 
Of course, if one chooses ip = i^i^), Eq-(12) can always be used to rewrite 
the relation in terms of p. For instance, in the case considered in Ref. [26] , ip 
proportional to nH^ implies that tfj scales with p^^~ . 



3 Inflation and Matter Creation 

As is well known, inflation is a theory of the early universe based not only 
on the possible existence of a primordial scalar fleld. To obtain the required 
dynamics, the potential of such a fleld needs to be able to, at least for a 
flnite period, generate a state of negative stress, which is the key condition to 
realize inflation. As shown by Guth and Sher|^l[], a prerequisite for inflation 



to work is a departure from thermodynamic equilibrium. In this context, 
it is naturally of interest to establish, at least qualitatively, how inflation 
can be implemented in the present irreversible matter creation theory. In 
order to analyze this issue, we flrst deflne an effective "adiabatic index" 
7* = (1 — 3^)7, so that Eq.(7) assumes the following FRW type form : 

RR+{^-^)R' + {^-^)k = , (13) 

which can be obtained using only the EFE (2) and (3) together with the 
effective equation of state p = (7* — l)p. As usual, for early times, we 
will neglect the spatial curvature contributions. In this case, Eq.(9) can be 
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rewritten as 



H+^H^ = Q . (14) 

Hence, from the expression for 7^, we see that the condition for exponential 
inflation [H = 0) is given by 

i) = 3nH , (15) 

or equivalently (from (8)), (3=1. Physically this is not a surprising fact, 
since the matter creation rate given by (15) has exactly the value that com- 
pensates for the dilution of particles due to expansion. As we shall see in 
section 5, exponential inflation occurs isothermically so that there is not 
extreme supercooling and violent subsequent reheating, as happens in all 
variants of inflation driven by a scalar field. As a matter of a fact, with con- 
tinuous matter creation, even in the power-law inflationary case, the decrease 
in temperature is much less than in the adiabatic case. The physical reason 
is quite simple. In this context, the entropy generation is concomitant with 
inflation, differently of what happens in the usual inflationary variants, where 
entropy is generated after inflation by a highly nonadiabatic process. In the 
new inflationary scenario, for instance, the temperature should decrease dur- 
ing the slow-rollover phase at least by a factor of 10~^^ in order to maintain 
the radiation entropy constant. In the present scenario, the matter creation 
continuously reheats the medium so that the temperature need not decrease 
so drastically. Essentially, this is the same result arising in the framework of 
inflationary models driven by bulk viscosity (see for instance Ref[|19|). 

Now, recalling that a violation of the strong energy condition (7^, < |) is 
a sufficient condition for "power law" inflation, one can show that (15), in 
this case, must be replaced by 

^>(l_A)3^i7 , (16) 
37 

or /3 > 1 — ^. Note that due to the matter creation process, either exponen- 
tial or "power law" inflation are now compatible with the existence of usual 
matter described by the "7-law" equation of state. In fact, Eqs. (15) and 
(16) do not impose any constraint on the 7 parameter. In this way, one may 
now refer, for instance, to a radiation or dust dominated "power law" infla- 



tion, as well as to different classes of de Sitter models (see also Ref. ||22|). In 
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addition, it is easy to see that the energy density, for each value of 7, scales 
with the temperature in the same fashion as happens for equilibrium states. 
For instance, in the era of radiation (7 = 4/3) one obtains p = aT"^. Only 
the time-dependence of the quantities is modified (see section 5, specially 
Eq.(42)). 

It should be observed that irreversible matter creation may also describe 
the so-called super (or pole) inflationary expansion [H > 0) in the terminol- 



ogy of Ref This kind of scenario appears, for instance, in string theory 
when taking into account the effects of the dilaton, a scalar field capable of 
driving superinfiation. Such a model has recently been proposed as a possible 



alternative to standard inflation (see Ref. ||3^ and references therein). In 



the present context, as one can see from (14), the condition H > will be 
satisfied if 7* < 0, that is, ip > 3nH. 

Another important point is related with the end of inflation, that is, the 
beginning of the FRW-type expansion. As the reader may conclude himself, 
the condition for inflation to come to an end can be obtained by the onset 
of violation of the above inequality (16). In fact, rewriting (14) as 

f = (l-^)i^^ , (17) 

it is clear that the stage of accelerated expansion will finish [R = 0) when 
ip = {1 — 2/3'y)3nH , that is, /? = 1 — 2/87. For instance, for 7 = 4/3 the 
universe evolves naturally from an exponential inflation to a FRW-type ex- 
pansion provided that the matter creation decreases in the interval 3nH/2 < 
ip < 3nH. For dust (7 = 1), a slightly different interval is required since the 
inflationary period will finish when ip = nH, that is, f3 = 1/3. Parentheti- 
cally, we remark that such conditions do not constrain the magnitude of the 
Hubble parameter to assume any specific value. As we shall see in section 6, 
this fact is closely related with the solution of the age problem which plagues 
the standard model for all values of the curvature parameter. 



4 The Simplest Class of Models 

Having in mind the choices previously made for ijj we now propose a specific 
matter creation scenario with a slightly modified creation rate. As remarked 
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in the introduction, we are not interested here in presenting a complete cos- 
mological scenario with matter creation, that is, a model describing the very 
early universe (including inflation) and the late stages of the evolution. Our 
goal is a much more limited one. We try to formulate a basic scenario or 
equivalently, a kind of hot big-bang model minimally modified due to the 
matter creation, over which the infiationary mechanism, as discussed in the 
later section, or any other process, may be further implemented. 

In our opinion, the simplest possible case, and probably the most phys- 
ically appealing too, at least for times later than the Planck era. is the one 
for which the characteristic time scale for matter creation is the Hubble time 
itself. Phenomenologically, this is equivalent (using (8)) to taking 



V' = SPnH , (18) 

where now /3 is a constant, which is presumably given by the particular 
physical model of matter creation. The above creation rate also simplifies 
considerably the task of solving eq. (13), since the effective "adiabatic in- 
dex" becomes 7* = 7(1 — (3) = const. In this case, it is readily seen that 
the generalized second-order FRW equation for R{t) given by (13) can be 
rewritten as 



RR + Ai?2 + AA; = , (19) 
the first integral of which is 

= . (20) 

where A = and ^4 is a positive constant (see eq.(2)). 

Using (20) one may express the energy density, the pressures (p and Pc) 
and the particle number density as functions solely of the scale factor R and 
of the P parameter. In fact, inserting (20) into (2), one obtains 

D 37(1-/3) 

P = Po(f) , (21) 

where po — SA/SnGR^* . The above equation shows that the densities of 
radiation and dust scale, respectively, as pr ~ and Pd ^ R ^ ' '. 

Hence, in a model with radiation and matter, the transition from radiation 
to a dust dominated phase, in the course of the expansion, happens exactly 
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as in the standard model. Note also that, although formally defined by the 
same expression (see (18)), the creation rates of radiation(7 = 4/3) and 
dust(7 = 1) are not equal, as they seem to be at first sight. For these cases, 
the above results are easily recovered in the usual manner, e.g. defining 
n = n.f + rifi and similar forms for the energy density and pressure. The 
corresponding dominant component then will determine the final form of all 
physical quantities (see also comment below Eq. (24)). 
Now, from (21), (5) and (6) 

D 37(1-/3) 

Pc = -hPo{^) , (22) 

with the total pressure Pt = p + Pc assuming the form 

/? 37(1-/3) 

Pt = {l*-l)p=b{l-f3)-l]Poi^) . (23) 

Finally, by integrating (4) with ip given by (18), or more directly, using 
eqs. (12) and (21), the expression for the particle number density can be 
written as 

/? 3(1-/3) 

n = n^Q) . (24) 

A clarifying comment about the meaning of equations (21)-(24) is now in 
order. Firstly, we note that (24) does not depend exphcitly on the 7 pa- 
rameter. Thus, the same scale law describes the evolution of the particle 
number density either for a dominant or a nondominant component. The 
effect of matter creation in both situations is measured by the (3 parameter. 
The situation is clearly different for the remaining equations, even though 
that those can also be applied for the nondominant component. In other 
words, for each expansion stage, the explicit time dependent form of the en- 
ergy density, equilibrium and creation pressures as well as the scale factor 
(sec Eq.(25)), depends exclusively on the dominant component, however, the 
creation of the other component is not completely supressed. Of course, this 
is the same kind of approximation commonly used in the standard FRW 
model. The only difference is that even in the radiation era, the dust compo- 
nent will have a nonvanishing creation pressure (see Eq.(22)), which although 
negligible in comparison with the creation pressure of radiation, will be re- 
sponsible by the baryon production in that phase. Such considerations will 
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be important when we discuss the thermodynamic behavior of these models 
(see section 5). 

As expected, for /3 = 0, eqs. (18)-(24) reduce to those of the standard 
FRW model for all values of the parameters k and 7 . In this case, the 
unified solution of (19) or, equivalently (20), was found by Assad and Lima 
[0 in terms of hypergeometric functions. Of course, such a solution can be 
adapted to the present case simply by replacing the "adiabatic index" 7 by 
the effective parameter 7*. 

For k = 0, the solution of (19) for all values of 7 and f3 can be written as 

2 

^^(1 — 0) 37(1-/3) 

R = Ro[l + ^^ ^ ^' Ho{t-to)] . (25) 

Note that in the limit 7^0 the above solution describes a de Sitter type 
universe for any value of j3. As remarked in the previous section, such a 
solution can now be obtained for /5 — > 1 and 7 assuming arbitrary values. 
Further, for flat models with 7^, > 0, we can choose to = 2i^~^/37(l — /?), so 
that (25) assumes a more familiar form, namely: 

^^(l ^ g) 37(1^-/3) 

R{t) = Ro r\ ^' Hot] . (26) 

If /3 = 0, (25) and (26) reduce to the well known expressions of the fiat FRW 
model. 

For k ^ 0, parametric solutions are usually more enlightening. By intro- 
ducing the conformal time coordinate, 

dt = Rdr] , (27) 

(19) can be recast as 

RR" + (A - 1)R'^ + MR^ = , (28) 

where the primes denote conformal time derivatives. 

Now, defining an auxiliary scale factor z = R'^, it is readily seen that 
(28) becomes 

." = ./ 7 = 5(J^ (29) 
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and 

z" + kAh = if , (30) 

whereas the first integral (20) is transformed into the energy conservation 
equation: 

1 '2 + 1^^2^2 ^ 1 2 2 (31) 

2 2 2 ^ ' 

where uj = and = Rl^. 

As one can see by direct substitution, the general solution of (30) or 
equivalently (31) is 

.^..-'"^1^1'" + ^ , (32) 

where 5 is an integration constant. Now, by choosing 5 = and using the 
inverse transformation R — z^^^, the general solution for the scale factor 
takes the following form: 

, (33) 



and 



Vk 



i(ry) = -R* y [ ^-/^^ -\ dr/ + const. (34) 

It should be remarked that the auxiliary scale factor z = R'^ shows the 
same dynamic behavior appearing in the standard FRW model, namely it 
evolves as a free particle {k — 0), a simple harmonic oscillator {k — 1) or 
an "anti-oscillator" {k — —1). In this sense, some basic characteristics of 
the standard FRW models arc not modified, namely open and flat universes 
expand forever whereas closed geometries exhibit a turning point cither when 
the universe expands away from the singularity {R = 0) or starts contracting 
from R — oo with the models presenting a big bounce. For completeness, we 
observe that the physical meaning of 2;* or equivalently i?* is readily obtained 
from the first integral (20). For instance, in the case of closed geometries with 
7* 7^ 2/3, R^ is just the value of R at the turning point, that is, for which 
i?(i?*) = 0. Accordingly, in the conformal time description we see from (31) 
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that the turning point 2;* (7* 7^ 2/3) corresponds to the amphtude of the 
related spring-mass system of unit mass (SHO). We leave it to the reader to 
verify that for singular flat models = Ro- 



5 Thermodynamic Behavior 

The matter creation formulation considered here is a clear consequence of 
nonequilibrium thermodynamics in the presence of gravitational fields [^2|, 
|23| , In this context, unlike other approaches to matter creation proposed 
in the literature (Cf. , for instance, |^ and [0), the explicit thermodynamic 



connection leads naturally to specific predictions on rates of variation of the 



entropy per particle and of the temperature. As shown in Ref. for the 
case of adiabatic matter creation, these rates are 

d = , (35) 

t = . (36) 

T dp „n 

where a is the dimensionless specific entropy and T is the temperature. 

For pedagogical convenience we first discuss the time dependence of the 
temperature. Using the 7— law equation of state, the temperature evolution 
equation (36) takes the form below. 



the integral of which is 



n{T) = n^i^y" . (38) 

o 

Now, replacing into (38) the deduced relation between p and n given by (12), 
the former can be written as 

P(r)=Po(^)'" • (39) 
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The above expressions for n{T) and p{T), are exactly the general expressions 



obeyed by a 7- fluid in the course of adiabatic expansion As remarked 
in section 2, for the case of radiation (7 = 4/3) the energy and particle 
number densities scale, respectively, as pr — and — so that the 
created radiation necessarily satisfies the usual equilibrium relations. This is 
a remarkable result. Differently from other approaches for matter creation, 
where expressions like (38) and (39) need to be assumed (see, for instance, 
Refs. 1^), the equihbrium relations here follow as a consequence of the 
"adiabatic" condition. In fact, as discussed in detail in the literature [^, ^], 
the condition (35) determines simultaneously the creation pressure form (5) 
and the temperature equation given by (36). On the other hand, combining 
p{R) given by (21), with (42), we obtain the following temperature law : 

p 3(7-l)(l-/3) 

T = U^) . (40) 

The above result shows us that exponential inflation {(3 = 1) occurs isother- 
mically regardless of the value of 7(see section 2). In fact, generically, the f3 
parameter works in the opposite sense of the expansion, that is, diminishing 
the cooling rate with respect to the case with no matter creation. In par- 
ticular, instead of the usual result, RT = const., valid for radiation(7 = |), 
we find TR^~f^ = const. Note also that by integrating (37) with n = the 
temperature law assumes a new form where the f3 parameter does not ap- 
pear explicitly, namely: N^^~"'^TR'^^^~^^ = const., which makes transparent 
the conclusion that if = const., the usual evolution law is recovered. This 
formula does not depend on the specific creation rate assumed in the present 
paper. In particular, for 7 = 4/3 one has 

N-^'^TR = const. , (41) 



as one should expect (see Ref. 0]). Now, recalling that for FRW geometries 
the frequency redshifts obeying, u ~ R~^, the above temperature law leads 
inevitably to the conclusion that the usual Planckian spectrum is destroyed in 
the course of the evolution, in particular, after decoupling. Nevertheless, as 
recently shown by one of us[0, for "adiabatic" matter creation the preserved 
spectral distribution is given by 
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. , (42) 

where iV(t) is the comoving time dependent number of photons and No is the 
constant value of evaluated at some fixed epoch, say, the present time. The 
above distribution is a consequence of the temperature evolution law as given 
by (41). When there is no creation, N{t) = No, and the usual Planckian form 
is recovered. In addition, it is readily seen that the equilibrium relations are 
recovered using such a spectrum. In fact, for 7 = |, it follows from (38) and 

(39) that n p^, and by introducing a new variable x = i^)^^, it is easy 
to see that 

POO 

p{T) = / pT{iy)du = aT* , (43) 
Jo 

where a is the usual radiation density constant. The spectrum given by (42) 
seems to be the most natural generalization of Planck's radiation formula in 
the presence of "adiabatic" photon production. More important still, (42) 
cannot be distinguished from the usual blackbody spectrum at the present 
epoch when we take T = To and N(to) = No- Therefore, models with "adi- 
abatic" photon creation may be compatible with the present isotropy and 
spectral distribution of the microwave background. Of course, since photons 
are injected satisfying (42), which is preserved in the course of the evolu- 
tion, there will be no distortions in the present relic radiation spectrum. 
Note that (42) is preserved precisely due to vahdity of the temperature law 
(41). In this concern, the macroscopic formulation adopted here seems to be 
naturally connected with some fundamental cosmological irreversible mech- 
anism (based on microphysics) , in which photons are quantum mechanically 
produced with the above thermal spectrum and baryons are asymmetrically 
created . 

Let us now consider the entropy behavior as defined in (35). Since 
cr = S/N, where S and are, respectively, the entropy of the dominant 
component and the corresponding number of particles, (35) can be rewritten 

as 

14 ■ («) 
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Hence, due to the matter creation processes, the universe does not expand 
adiabatically as happens in the standard FRW models[^]. Besides, since up 
to a constant factor one has = nR^, it follows from (24) that increases 
as a power of R, that is. 



Further, from eq. (37), S = So{N/No), and using the above expression 
one may write the photon entropy as (from now on indexes r and b refer, 
respectively, to radiation and baryon component (dust)) 



where 5*0,. ~ 10*^^ is the present observed radiation entropy (dimensionless). 
As happens in the standard model, although remaining nearly constant, the 
specific entropy defined by (35) does not play an important physical role. 
As we know, some physically meaningfuU informations, as for instance, in 
nucleosynthesis studies as well as for the structure formation problem, are 
encoded in the specific radiation entropy per baryon. Such a quantity, de- 
fined by arb = is proportional to the photon-baryon ratio and, up to 
short aniquillation period, also remains constant in the standard model. In 
the present context, since photons are thermally produced and baryons are 
continuously created, a nearly constant behavior of arb should also be ex- 
pected. In fact, the net number of baryons in the comoving volume is given 
by (see discussion below Eq.(24)) 



where Nob is the present baryon number. It thus follows from (46) that 
arb = tF' = '^o ~ 10^°, is the present photon to baryon ratio. As a kind of 

ob 

consistency check, we notice that if one write the specific entropy(per baryon) 
in the usual form 



N 




(45) 




(46) 




(47) 



arb — 



4aT^ 
Srib 



(48) 
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the above result is recovered since Tj} and Ub evolve following the same scale 
law(see Eqs.(40) and (24)). Therefore, the adiabatic formulation does not 
provide any explanation for the present value of photon to baryon ratio. As 
happens in the big-bang model, the value of cTq is just an initial condition. In 
the present model, particles are created in spacetime with the same temper- 
ature as the already existing ones have, otherwise the specific entropy could 
not remain constant, as happens in the more general formulation proposed 
in Ref [^. Naturally, from a thermodynamic point of view the model is ir- 



reversible. The burst of entropy is closely related with the creation of matter 
and radiation. 

We would like to stress that all important thermodynamic results of the 
standard FRW models like S = So, N = No and the radiation temperature 
scaling T cx are recovered for /? = 0. Finally, it is interesting to remark 
that the models presented here may significantly alter the standard predic- 
tions of cosmic abundances, since they alter the expansion rate and predict 
a new temperature law. Such results points to possible limitations on the f3 
parameter imposed through constraints from nucleosynthesis. This issue will 
be addressed elsewhere. 



6 Some Observational Aspects 

Now we illustrate some observable predictions of the models proposed in the 
preceding sections. Following standard lines we define the physical parame- 
ters q = (deacceleration parameter), H = ^ (Hubble parameter) and 

f2 = ^ (density parameter), where pc = |^ is the critical density. 
Inserting the above quantities into Eqs. (2) and (19) we have 

and 

A = (fi_i)if2 _ (50) 

Therefore, it is clear that if > 1, that is, if g > the universe is 

positively curved with p > pc, whereas if g < ^"^(-^"^^'^ it is negatively curved 
or flat, respectively, with p < pc. For = the usual expression for FRW 
models are recovered. However, the positivity of Vt does not restrict q to be 
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positive if 7 > 2/3 as happens in the standard universe. To be more specific, 
in the present dust phase (7 = 1), the above expressions reduce to 

Q^^^^ , (51) 
^ -(fio-l)/// . (52) 



2 



Ro 

Therefore, if Qo = ^-^^ we have = 1 and from (52), the presently observed 
universe is flat. However, regardless of the value of k, the deacceleration 
parameter may be negative since the constraint fio > can be satisfied for 
5o < 0, provided that /3 > 1/3. As we shall see next, such a fact allows us to 
solve the age problem in the present context. 

As we know, the age of the universe is found by integrating the generalized 
first integral (20). By expressing the constant A in terms of fig, Ro and 
it is straightforward to show that 

(^) =H',[i-n^ + nA ] , (53) 

the solution of which may be expressed as a formula for the time t in terms 
of R, 

fR/Ro 

t-U = Ho'^ [l-Qo + ^oX-'^^Y^'^dx , (54) 

where is the time for which R = R^. 

For singular models, the present age of the universe is defined by taking 
R^ = t{R^) = so that it is given by 

t = Ho'^ [\l-no + ^oX~^^]~^^^dx , (55) 

Jo 

which for P = 0, that is, A = is exactly the same as in the standard 

model. In Figs. 1 and 2 we show the age of the universe in units of for 
the cases of dust and radiation dominated universes and some selected values 
of (3. Observe that the universe can be old enough even when considering a 
radiation dominated phase today, as has been suggested sometimes. Another 
important conclusion is that if /3 > 1 — 2/37 and < < 1 the oldest 
universe is the flat one {flo = !)• this context, our matter creation ansatz 
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(18) changes the predictions of standard cosmology in such a way that it 
solves the problem of reconciling observations with the inflationary scenario. 
Note also that open models with a small parameter are also ruled out by 
recent data regardless of the value of 13. 

The solution of the age problem and its noticeable dependence on the (5 
parameter can be exactly determined for fiat models. In this case (fio = 1), 
one can see that the age parameter reduces to 

*»'»=3Ml^ ■ '''' 

In a matter-dominated universe we have then Hoto = 2/3(1 — (3), and in this 
form it is easiest to see how matter creation could solve the age problem 
suggested by the latest direct measurements of the Hubble constant done 
by Pierce et and Freedman et who found, respectively, Ho = 

87 ± 7Kms~^Mpc~^ and if^ = 80 ± llkms"^ Mpc'^ . Assuming no matter 
creation (/? = 0), these values of Ho imply that the expansion age of a dust- 
filled, fiat universe would be about either 7.3 x 10^ years or 8.2 x 10^ years, 
in direct contrast to the measured ages of some stars and stellar systems, 
believed to be at least some (16 ± 3) x 10^ years old or even older if one adds 
a realistic incubation time |4^. Such measurements restrict the parameters 
Hoto to the following intervals (P and F reffer to the values of Ho given, 
respectively, in Refs.pH] and p9|l) 



1.09 < Hoto < 1.86 (P) (57) 

and 



0.85 < Hoto < 1.91 (F), (58) 

when in the standard fiat model (7 = l,/5 = 0), one would obtain exactly 
2/3. In this fashion, these recent measurements point to a serious crisis of 
the standard model. It would not be surprising if expectations that it could, 
at least marginally, be compatible with the age of the universe, will gradually 
be forgotten (see Figs.l and 2). 

As can be easily seen from (56), matter creation naturally solves this 
problem, increasing the parameter Hoto while preserving the overall FRW 
evolution scheme. From (56), (57) and (58), the constraints on the 13 param- 
eter for both cases are readily computed to be 
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0.38 <p< 0.64 (P) (59) 

and 

0.21 </3< 0.64 (F). (60) 

Note that due to the error bar in the values of Ho, the upper bound of /J 
does not depends on the particular set of measurements. Naturally, a more 
realistic range of (3 probably will have an upper bound slightly smaller than 
0.64 and a lower bound between 0.21 and 0.38. This will happen, for instance, 
if the future measurements of Ho improve by at least one order of magnitude. 
Parenthetically, there are some independent indications requiring the present 
expansion rate to lie in the range Ho = 80 ± 5Kms~'^Mpc~^ (see and 



references there in). With such a precision, the P parameter will fall on the 
interval 0.34 < /3 < 0.60. 

As remarked earlier, the solution of the age problem is due to the pos- 
sibility of go < nowadays. For instance, from (49) with fio = 1 we can 
rewrite the result (59) as Hoto = 1/(1 + go)- For P in the above interval, 
we have for a dust filled universe —0.44 < go < —0.1. In connection to this 
we note that cosmological constant models with a reasonable value of fl\ 
are also believed to solve easily the "age problem". However, as recently 
shown by Maoz and Rix the computed rate of gravitational lensing in 



such models constrain severely f^A when confronted with the existing lensing 
observations. Analogously, such a result seems to point out to similar limits 
on the /3 parameter, thereby leading to values of Hoto below the range given 
by (57) or (58). In this concern, we remark that models with matter creation 
behave like scenarios driven by a decaying Aterm, instead of models with 
cosmological constant [^, Moreover, for this kind of models, a lower 
lensing rate is usually predicted since the distance to an object with redshift 
z tends to be smaller than the distance to the same object for a model with 



A constant [^| 



As we have discussed in the previous section, photons are always created 
in equilibrium with the already existing radiation. However, as we have seen, 
they are not dominant nowadays. In the present dust phase, the matter 
creation rate is given by 

iJo = ^PnoHo . (61) 
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For Uo ~ 10^® nucleons cm~^ and ~ 10^°|/r we have ipo ~ 10"^^ 
nucleons cm~^yr~^, which is nearly the same rate predicted by the steady- 
state universe ^ and also by some decaying A cosmologies PT|. Of course, 
this matter creation rate is presently far below detectable limits. 



7 Conclusion 

Almost all research efforts related to the physical processes in the cosmolog- 
ical domain have been devoted to the standard model in its different phases. 
As we know, such a model is thermodynamically characterized by two differ- 
ent although related features, namely: entropy conservation (5." = 0) and 
number particle conservation {N!^ = 0), where S°' and A^" are, respectively, 
the four-vectors of entropy and number of particles. 

In this paper we have investigated some cosmological consequences arising 
when one changes the second (and consequently the first) of the above prop- 
erties. Of course this is not new as may easily be observed in the extensive 
literature on this subject (see p| and references therein). The new fact jus- 
tifying the present work is that we have considered a recent thermodynamic 
approach for which matter creation, at the expense of the gravitational field, 
has been properly constrained by the usual requirements of nonequilibrium 
thermodynamics [22-25]. In this context, thermodynamic predictions such 
as the temperature law and the variation rate of the entropy are computed 
from first principles. Such an approach points to a possible revival of interest 
in models with radiation and/or matter creation, which is one of the most 
challenging problems of theoretical physics. 

For "adiabatic" matter creation a general expression relating the parti- 
cle number and energy densities which holds regardless of the creation rate 
has been deduced and the conditions for generating inflation have also been 
established. It was also shown that a minimally modified big bang model 
with adiabatic matter creation can be easily implemented and predicts in- 
teresting cosmological consequences. In fact, for a creation rate = SjSnH, 
the known perfect fluid solutions can be adapted by using a slightly modified 
"adiabatic index" 7,, = 7(1 — /?). The thermodynamic behavior is also readily 
computed making clear how the above mentioned properties of the standard 
model are quantitatively changed. In particular, instead of constant entropy 
and T ~ R-\ we found T ~ R-^^-^^ and S ~ R^^. In addition, since the 
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specific entropy of radiation is very liigli, one can sliow tliat as long as matter 
is in tliermal contact witli radiation it will follow the same temperature law 
as the radiation, so that the thermal history is accessible as in the standard 
model. The model is also able to harmonize a FRW-type picture with the 
discrepancy existing between the latest measurements of the Hubble param- 
eter and the age of the universe, as predicted by the standard model. Of 
course, in order to have a viable alternative to the standard FRW model, 
other well known cosmological tests need to be investigated. Further details 
of our model will be published elsewhere. 
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Captions for Figures 



Fig. 1 



Fig. 1 - The age parameter of a matter dominated universe as a 
function of Clg for some selected values of /?. The horizontal dotted 
and sohd hnes indicate the latest observational constraints (see Eqs. 
(52) and (53)). The solid curves represent two limiting cases, the 
standard model (/? = 0) and a de Sitter type universe (/? = !). Curves 
A, B and C respectively represent models with (3 = 0.44, (3 = 0.54 
and P = 0.63. According to Pierce et al. data, the standard dust 
model is ruled out regardless of the value of JIq whereas for Friedman 
et al. only extremely open FRW models may be compatible with the 
observations. For /? > 1/3 the matter creation process rehabilitates 
the fiat universe, as predicted by inflation. 



Fig. 2 



Fig. 2 - The same graph of Fig. 1 for a radiation dominated universe. 
The solid curves now indicate the radiation filled FRW model (/3 = 0) 
and a de Sitter-like radiation universe {(3 = 1). Generically, the curves 
are displaced downwards in comparison with the case of dust so that 
the models are compatible with data only for higher values of (3. 
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